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Abstract 

The Kroemer model for the Gunn effect in a circular geometry (Corbino disks) 
has been numerically solved. The results have been interpreted by means of 
asymptotic calculations. Above a certain onset dc voltage bias, axisymmetric 
pulses of the electric field are periodically shed by an inner circular cathode. 
These pulses decay as they move towards the outer anode, which they may not 
reach. As a pulse advances, the external current increases continuously until 
a new pulse is generated. Then the current abruptly decreases, in agreement 
with existing experimental results. Depending on the bias, more complex 
patterns with multiple pulse shedding are possible. 
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I. INTRODUCTION 

Propagation of pulses naturally occur in excitable media, that exhibit a large response 
when a sufficiently strong perturbation disturbs the only stable stationary homogeneous 
state UJ. Examples are the propagation of an action potential along the axon of a nerve 
0, the propagation of a grass fire on a prairie, pulse propagation through cardiac cells 
0, reaction-diffusion || or ecological systems |T[]. While a vast literature is devoted to 
the mathematical description of pulses propagating in unbounded media, less is known 
about pulse generation from boundaries and propagation in finite domains, particularly in 
multidimensional domains. 

That boundaries and boundary condition play an important role in pulse creation and 
annihilation is well understood in the subject of semiconductor instabilities [|J. Device 
geometry and bias conditions are crucial for instabilities and related nonlinear dynamics to 
appear. Nice examples can be found in the experiments by Willing and Maan on repeated 
pulse propagation in semi- insulating GaAs f5[. They considered rectangular samples with 
two attached parallel planar contacts or with point contacts at different dc voltage. In the 
first situation, planar pulses were periodically generated at the cathode and moved towards 
the anode where they disappeared. This phenomenon is analogous to the well-known Gunn 
effect in bulk n-GaAs [|J. In the case of point contacts, circular waves were repeatedly 
generated at the cathode and vanished before arriving at the anode ||. Theoretical studies 
of self- sustained oscillations in semi-insulating GaAs are scarce even in one- dimensional 
geometries (cf. Ref. J/J and references cited therein). However, the observed phenomena 
can be qualitatively understood within the simpler Kroemer model for the Gunn effect in 
bulk n-GaAs |§ . An asymptotic study of this model on a one-dimensional spatial support 
can be used to understand pulse propagation in samples with planar contacts |J. A simple 
study of pulse propagation in samples with point contacts could consist of analyzing the 
Kroemer model in a axisymmetric sample: a circular sample of bulk n-GaAs with a point 
contact (cathode) at its center and an attached concentric circular outer contact (anode). 
This configuration is known as Corbino geometry In this paper, we carry out a numerical 
study of pulse dynamics in the Kroemer model with Corbino geometry. Depending on the dc 
voltage bias and contact resistivity, we observe stationary field and current and self-sustained 
oscillations (periodic or not) due to pulse propagation and recycling. Pulses may or may 
not arrive at the anode before a new pulse is generated at the cathode. These results are 
presented in Section |TTT| after a short description of the Kroemer model in Corbino geometry 
is given in Section |I[ The numerical results are interpreted by means of an asymptotic 
analysis in Section [IV]. Section |V| contains our conclusions and the numerical method we 
use is described in Appendix [A]. 



II. EQUATIONS AND BOUNDARY CONDITIONS 

The Kroemer model consists of the following equations and boundary conditions (in 
dimensionless units) for the concentration of free carriers (electrons), n, and the electric 
potential, ip: 
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Here ([!]) and (||) are the charge continuity and Poisson equations, respectively. The dimen- 
sionless electric field is E = Vy? and E = \E\. In these equations, the electron density 
has been scaled with the uniform concentration of donor impurities in the semiconductor, 
No = 10 15 cm -3 , and the electric field with the field characterizing the intervalley transfer 
responsible for the negative differential mobility involved in the Gunn oscillation, Er = 3.1 
kV/cm. Distances and times have been measured with the dielectric length and the dielec- 
tric relaxation time, l\ = eE R /(eN D ) « 0.276/im, li/(p Q E R ) « 1.02 ps, respectively (/i is 
the zero-field electron mobility; see, e.g., |10| for details). The unit of electric potential is 
ErIi ~ 0.011 V. The carrier drift velocity of Eq. (^|), v(E), is already written in dimension- 
less units, and it has been depicted in Fig. |l|. We assume that the diffusion coefficient is 
constant, S « 0.013 (at 20K). 

Boundary and bias conditions need to be imposed at the interfaces between semiconduc- 
tor and contacts, S c a . Our boundary conditions @ and (f|) assume that the normal com- 
ponents of electron current density and electric field are proportional at the semiconductor- 
contact boundary (Ohm's law) JL0|, (in these equations, N is the unit normal to S c a , directed 



towards the semiconductor). For simplicity, we choose all contact resistivities p to be equal. 
Bias conditions are chosen to be (p — at the cathode S c (injecting contact) and <p = $ 
(the applied voltage) at the anode S a (receiving contact). Typically 5 > is very small, 
so that diffusion matters only inside boundary layers near the contacts or inside thin shock 
10, 9| . The latter are charge accumulations that will be treated simply as disconti- 
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Kj. Thus diffusion effects may be left out of the conservation 

0, the first boundary condition in 



nuities of the electric field 

equation (|l|) when interpreting the results. If we set 5 
(0) should be omitted. 

We can write an Ampere's equation for the total current density (electronic plus dis- 
placement), j, by eliminating n from (0) using (0): 
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In the Corbino geometry considered in this paper, this equation can be simplified further. 
Let r c and r a > r c be the radii of cathode and anode, respectively. Electric field and current 
density are now directed along the radial direction, E = E(r,t)r/r, E(r,t) = dip{r,t)/dr 
and j = J(t)r/r 2 , so that (|6]) becomes 
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where 2nJ(t) is the current through the external circuit, i(t) = J* E j ■ N dA = 2nJ(t). 
Equation (|7|) for E(r,t) and J(t) should be solved with the following bias and boundary 
conditions 
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(8) 

E = P \z ~ ~ZT ) at r = r C} r a , (9) 

where L = r a — r c and = 

It is known (see, e.g., fllP"! ) that planar dipole waves may appear in long samples when 
p > 4/3, for which the straight line j = E/p representing Ohm's law intersects the drift 
velocity curve j = v(E) at a point (E c ,j c ) on the second branch of this curve, as sketched in 
Fig. |l| (In particular E c rs p 1//4 and j c ~ p -3 / 4 for large p). This is the case we will consider 
in the present paper. 



III. NUMERICAL RESULTS 

We have solved numerically Equations (^) and (§) together with the boundary conditions 
([]) at r = r c and r = r a . Appropriate initial conditions were given for E(r, 0). Parameter 
values were p = 2, 5 = 0.013, r c = 10 and r a = 50 and 90 (i.e. L = 40 and 80 resp.). The 
bias was used as a control parameter for two different electron velocity curves v s = and 
v s = 0.1, representing zero and nonzero saturation velocities at high electric fields. These 



parameter values are appropriate for n-GaAs and consistent with previous studies |L0 
Consider firstly the characteristic current-voltage curve J(<f>) of Fig. |2|. There we can mark 
three different regimes, already present in experiments 0. 

Regime I: < <fi < <p a ~ 0.168. Stable solutions are stationary and J(4>) is well 
approximated by a straight line with slope ~ 31.5. 

Regime II: <p a < < 30 a /2 « 0.25. Above the onset bias for current oscillations, there 
are small-amplitude (10-20% of the overall current signal) sinusoidal current self-oscillations. 
The oscillation maxima and minima are about J c ~ 5.4 and J m in ~ 4, respectively. The 
electric field profile is a triangular pulse which is recycled at the cathode, advances and it 
soon disappears at r w 25 (quenched-mode oscillation). See Fig. || 

Regime III: 30 a /2 < < W . The upper critical bias W is finite for v s > (e.g. 
0^ ~ 15 if v s = 0.1) and infinite for v s = 0. There are large-amplitude (60% of the 
overall current signal) current self-oscillations. These oscillations are mostly time-periodic, 
although there are narrow bias intervals of aperiodic oscillations. Their maxima J m ax(4>) 
are always close to J c , while their minima J m in{4>) take values on a wider range of currents. 
Electric field profiles consist of moving triangular pulses. For > 0^, the stable field 
profile is again stationary as shown in Fig. [|. These results qualitatively agree with the 
experimental observations of self-sustained oscillations in semi-insulating GaAs reported in 
Ref. ||. Experimental samples were rectangular and contained two well-separated point 
contacts. It was observed that self-oscillations of the current were due to circular dipole 
waves that were recycled at the cathode, expanded towards the anode and vanished without 
ever reaching it. The current signal was similar to that in Fig. |5| below. Semi-insulating GaAs 
is described by model equations different from the Kroemer model. However, qualitative 
agreement of experimental data with our results for the axisymmetric Kroemer model suggest 



that a theoretical interpretation of self-oscillations similar to that in Section [TV] could also 
be appropriate for semi-insulating GaAs. 
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We shall now describe the more interesting self-oscillations in Regime III, starting with a 
bias interval of time-periodic oscillations. Stationary solutions and their stability properties 
will be described elsewhere. 



A. Time-periodic oscillations for 0.35 < (ft < 0.5 and v s = 

Fig. |5] shows one period of J(t) for different bias values in this interval. The electric field 
profile E(r,t) consists of a single triangular pulse traveling towards the anode when J(t) is 
increasing. When J(t) decreases, one triangular pulse disappears and a new one appears at 
the cathode. The gradual increase of the current when there is only one pulse in the sample 
lasts longer than the drop to low current values. Notice that the current drop lasts the same 
for all bias values, while the stage of current growth increases with (ft. 

This situation is remarkably different from that for the one-dimensional geometry cor- 
responding to parallel planar contacts: the current signal is flat when there is only a single 



pulse far from the contacts and the stages of current increase and drop are very short [JTO 
Other noticeable features in Fig. [S] are: (i) a new wave is nucleated as J surpasses a critical 
value J c (bias independent); (ii) the current overshoot above J c decreases as (ft increases and 
(iii) there is a second local maximum of the current and the width of the region between 
the two local maxima increases with (ft. Let us explain in more detail the field profiles 
corresponding to these stages of the self-oscillation. 

Fig. |6| and Fig. |7| show details of the current signal during one period of the self-oscillation 
for (ft = 0.38 and v s — 0. Also shown are the maxima of the electric field the different pulses 
that appear in this time interval. From t — to t — 195, there is a single pulse moving 
towards r a . This wave is roughly a straight triangle of height and width E + (t), which is the 
maximum of the field inside the wave at time t. The back of the pulse can be approximated 
by a shock wave located at Rbif). As we describe below, self-oscillations in the bias interval 
considered here involves creation of new pulses and of transient field disturbances at the 
cathode. In Fig. [7], R' b (t) and Rf %s {t) are the locations of the maximum field of a new pulse 
and of a transient pulse-like field disturbance (both shed at the cathode), respectively. Fig. 
||(a) shows the electric field profile at the times marked (1) and (2) in Fig. |. Fig. ||(a) 
depicts the electric field near the cathode. A new pulse is shed from the cathode when 
J surpasses the critical value J c at time (4) in Fig. || Between times (3) and (4) (with 
J < J c ), a field disturbance is shed from the cathode and it shrinks rapidly as it advances. 
Eventually as the first pulse disappears (slightly after the time corresponding to the second 
local maximum of the current), the new pulse shed at time (4) becomes the only one in the 
sample. A new oscillation period starts then. Fig. |S] shows details of the field profiles when 
there is more than one pulse in the sample. 



B. Aperiodic current self-oscillations 

For biases (ft > 0.5, periodic self-oscillations of the current alternate with voltage ranges 
of aperiodic oscillations. Fig. [T^ shows that the corresponding current signals may be rather 
complex, with several maxima and current overshoots above J c . 
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Typically several pulses are present during different time intervals of the self-oscillation. 
The leftmost pulse may shrink as the second pulse increases, or it may reach it and coalesce 
with that pulse. Meanwhile the rightmost pulse may shrink or reach the anode. Pulses may 
be shed from the cathode or nucleate inside the sample. Resulting current signals may even 
be apparently chaotic. A detailed study of all cases that are possible depending on the bias 
will not be attempted here. Fig. |TT] depicts the local maxima of the current as a function of 
bias. Loss of periodicity at narrow bias intervals is apparent. 

Comparing the case v s > to that with v s = 0, we observe that the pulses move faster, 
the oscillations have smaller amplitudes and it is easier for several pulses to coexist when 
v s > 0. This may result in more complex shapes of the current signal, as shown in Fig. [12. 
The current signal is periodic in Figures [121(a) and |l2|(b) although the period is longer in 
the latter due to separated current bursts. The current signal in Fig. |12|(c) is aperiodic. 



IV. INTERPRETATION OF THE NUMERICAL RESULTS 

In order to understand the shape of J(t) in the time periodic regime and the fact that 
the pulses may vanish before reaching the anode, we present in this section a qualitative 
description of the asymptotic oscillatory solution for r a 3> r c 1. In this case there is a 
wide range of voltages for which the pulses are detached from the contacts during most of 
their evolution time. A full analysis of the pulse dynamics should include descriptions of the 
evolution of a pulse far from the contacts and of the generation of new pulses at the injecting 



contact. The latter process is essentially as for the planar case |T0| when r c ^> 1, because the 
effect of the geometrical divergence is then negligible around the cathode. Analysis shows 
that, roughly, a pulse is shed when the current increases to J = J c « j c r c , and then the 
current decreases while the new pulse grows and separates from the cathode. These results 
are in line with the numerical simulations of the previous section. 

In the remainder of this section we focus on the evolution of a pulse detached from the 
cathode. As in the planar case, such a pulse is a straight triangle made of a trailing edge 



which is a shock and a leading ramp which is a region depleted of electrons ]10 . 

Consider first the region outside of the pulse. Time and space derivatives can be neglected 
in (0) for this region, which covers most of the sample, leading to the appoximate solution 
v (E) = J/r, which implies E(r,t) = Ei(J/r), where E\(j) < -E^C?) are the two solutions of 
v(E) = j for v s < j < %. If J <C r, the first branch of v(E) is linear, and we have E(r, t) w 
J/r. The area under this stationary field profile is § out = f£ a E^J/r) dr « J\n(r a /r c ). 

The speed of trailing edge of the pulse, at r = Rb(t), is given by the equal area rule for a 
shock rising the field from £L to E + ; i.e. V(E + , EJ) = v(E) dE/(E + — EJ). For large 
pulses, E + ^> 1 and the field immediately at the left of the shock is E- < Em = 0(1), so 
that we may approximate 

V(E +1 E^) = v s + E+ _ E Vs+ E^ ' 

with 

/•CO 

C= / [v(E)-v s ]dE. 

JEi(vs) 
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Here we have used that V ~ v s and EL ~ Ei(v a ) as E + 3> 1. If u s = 0, then C = ir/4 and 
V(E + , EJ) ~ 7r/ (4E + ). In this case, the trailing front velocity is small and small waves move 
faster than large ones. If v s > 0, then the waves move at a speed close to the saturation 
speed v s . 

The electron density at the leading ramp of the pulse is almost zero, so that the field 
obeys the Poisson equation 1 + r _1 d(rE)/dr — n ~ 0, whose solution is 



r 



2 



E(r,t)= 2^ — ~ r -- r (n) 

in the pulse. Here the constant of integration r w (t) is the intersection between the prolon- 
gation of the ramp and the r axis, and use has been made of the condition that the width 
of the pulse W = r w — Rb ~ E + satisfies r w ^> W ^> 1 to simplify the result. The area of 
the triangular pulse is $j n ~ E^_/2, and therefore the bias condition (|8|) becomes 

$ = $ in + $ out ~ & + J In (^j . (12) 

The speed of the ramp, dr w /dt, coincides with the speed of the electrons immediately 
ahead of the pulse 

dr w J fio\ 

~di ~T W [L6) 

This result can also be obtained by inserting (|TT|) and n = into (^). Equations ( p!2[) and 
(PH), along with 

for the speed of the trailing shock (at Rb ~ r w — E + ), suffice to determine the time evolution 
of J, r w and E + from a given initial state. Equation (0) can be rewritten as 

dJ 1 W^-,Jx 

(15) 



ln(r /r c ) I Vr 



$ - Jin 

where use has been made of (|1^) and (|T3|). 

The problem can be further simplified noticing that the speeds of the shock and of the 
ramp are nearly equal to each other during most of the pulse lifetime, 

". + §- = f (16) 

This occurs because the width of the pulse is small compared with the total distance it travels 
in the sample. Then a mismatch of the two speeds would lead either to the disappearance 
of the pulse or to its growth above the maximum size allowed by the bias in a time of 
order E + / max(V, J/r w ). This time is short compared with the pulse lifetime. An algebraic 
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relation between J and r w can be obtained eliminating E + between (|T|) and (|T^). This is 
plotted as a dashed curve in Fig. [l^, which is the phase plane of ( |13|) and ( |l5l) for a case 
with v s = 0. As can be seen, the trajectories tend rapidly to the lower branch of the dashed 
curve and then rise along this branch obeying ( |13"D until either: (i) r w reaches its maximum 
possible value, corresponding to the turning point T, or (ii) J reaches the critical value J c 
for nucleation of a new pulse, whichever happens first. 
When v s = the turning point is 



Ji 



2$ 



31n te) 71-111 fe) 

Then maximum radius of the pulse is either r w j- or 



2$\ 2 



(17) 
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if J c < Jt- The time dependence of J from its minimum value during an oscillation period 




can be obtained by integrating ( |T3"D along the dashed line (nullcline) of Fig. |H| In the 
present case of v s = 0, Eq. (|13|) can be written as [$ — (3 J/2) In (r a /r c )] dJ/dt = 7r 2 /32, with 
the help of ( |T2| ) and (p!6|). Upon integrating this equation, 

J(t) = J T fl- A fl-^L) -±| , (19) 

(20) 

We have used as initial condition the minimum value of the current during one oscillation 
period, J(0) = J m in- A comparison of this approximation to a direct numerical solution of 
the whole problem is shown in Fig. |TJ|. 

For shorter samples such that r a < r wc and large bias, the pulse reaches the anode before 
its maximal radius. Then the situation resembles the one- dimensional Gunn effect with large 
voltage: a pulse reaches the anode before a new pulse is nucleated at the cathode |10| . If we 
relax the assumption v s = 0, we find that the pulse speed is faster according to (|T0|) . As in 
the case of v s = 0, there is also a maximum radius of the pulse (for large enough samples). 
In fact, the nullcline dJ/dr w = has a turning point which can be found by solving 



s ~3 



for z > as a function of $. Then Jt and t w t are given by 



z = E 4 



$ — J T In ( — 



Jt z 



(21) 

(22) 
(23) 



S 



An approximate solution of the system (|T3| ) and flT5|) can be found as before. The maximum 
radius can be found by solving the simplified problem with initial condition J min and calcu- 
lating the time that J needs to reach J c . In general, the pulse reaches either its maximum 
radius or the anode (depending on L and 0) much earlier than in the case of v s = 0. This 
leads to small-amplitude current oscillations which may be rather irregular because J may 
be above J c quite often, thereby producing new pulses near the cathode. See Fig. [12. 



V. CONCLUSIONS 



We have studied numerically repeated generation and motion of axisymmetric pulses in 
a two-dimensional n-GaAs sample with a Corbino geometry. The field inside these pulses 
decreases as they advance and expand, so as to compensate their larger extension. Simulta- 
neously, the current increases until a critical value is reached and a new pulse is triggered at 
the central point contact. The current signal presents different patterns depending on the 
applied dc voltage bias. Just above the onset for self-oscillations, their amplitude is small 
and the pulse dies off shortly after it is generated. For larger voltages, self-oscillations have 
larger amplitude and pulses may or may not reach the outer sample boundary depending 
on the size thereof and bias. For sufficiently large samples, the pulse radius cannot surpass 
a maximum value given by an approximate analytical formula. Regions of aperiodic oscilla- 
tions due to multi-pulse dynamics are interspersed with more regular periodic oscillations. 
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APPENDIX A: OUTLINE OF THE NUMERICAL METHOD 



We have used an efficient numerical scheme for partial differential equations with an 



integral constraint described and proved to converge in Ref. Radial derivatives are 

approximated by central differences and a first-order implicit Euler method is used to inte- 
grate the resulting differential equations in time. This procedure results in having to solve 
a system of (N + 1) linear equations for the values of the electric field and J at time t n+ i 
in terms of their previous values. The block matrix formulation of this system is shown in 
Fig. |T^. There T, u and v are a N x N tridiagonal matrix, a 1 x N row vector and aiVx 1 
column vector, respectively. Our system is therefore equivalent to 

T • E + Jv = s, (Al) 
u • E + Ja = 0. (A2) 

This system can be efficiently solved by solving the following two systems with the same 
tridiagonal matrix: 

T ■ y = s, (A3) 
T z = v. (A4) 
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In terms of y and z, we obtain 



E 



Jz, 



J 



uy 



u • z — a 



(AS) 
(A6) 



Thus we proceed by firstly obtaining the LU factorization of T and secondly carrying out 
two back substitution processes to solve (|A3|) and (|A4[). Then (|AB|) and ( |A5| ) yield J and 
E, respectively. 
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FIGURES 



Ohm's law 




FIG. 1. Drift velocities and Ohm's law. v(E) = \v(E)\ has a maximum vm = 

3 3/4 /4 

at 

E = E M = 1/3 1/4 (for v s = 0), followed by a region of negative differential mobility for E > Em- 
At large fields the electron velocity monotonically decreases to a value v s , which may be 

zero. 
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FIG. 2. Current-voltage characteristic curve J{4>) for v s = 0. If < <fr < (f> a , the sta- 
ble electric field profile is stationary. For <fi > cf) a , we have depicted the maxima, minima 
and time-averaged value of the current self-oscillations. There are small-amplitude current 
self-oscillations if (j) a < <f> < 2cf) a /3, and large- amplitude self-oscillations due to pulse recycling 
and motion if 2</> Q /3 < <p < 4> w ). 4> w is finite for v s > and infinite for v s = 0. 
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FIG. 3. Total current density (left) and electric field profiles (right) for <j> = 0.18. The electric 
field profiles are depicted at the times marked on the graph of J(t). The horizontal line in the 
latter corresponds to the value J c . 
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FIG. 4. Stationary solution for <fr = 20, L = 80 and v s = 0.1. J evolves towards J ~ 6.2 greater 
than J c = 5.77. The maximum electric field, E + « 26.75, is reached at Rf> 35.12. Field values 
at the boundaries are E(r c ) = 1.25 and E(r a ) = 0.14. 
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FIG. 6. Evolution of J(t) and E + (t) during one oscillation period for cf) = 0.38 and v s = 0, 
during the time interval < t < At ss 200, where there is a single pulse in the sample. Marked 
times are: (1) 30, (2) 180, (3) 195, (4) 215. 
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FIG. 7. Evolution of J(t) and E + {t) during one oscillation period for (ft = 0.38 and v s = 0. (a) 
Stage in which there are multiple pulses, (b) Rb(t) during the multi-pulse stage. R^ iS (t) and R' b (t) 
are the locations of the maximum fields of the transient field disturbance and of the new pulse, 
respectively. Marked times are: (3) 195, (4) 215, (5) 230, (6) 248, and (7) 253. 
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FIG. 8. (a) Electric field profile during the the single-pulse stage at times (1) and (2) of Fig. 
^. (b) Details of the field profile near the cathode: the slope at r c increases with J; dE{r c )/dr = 
for J = J c k, 5.24 at time (3). 
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FIG. 9. (a) Electric field profiles during the multi-pulse stage at times marked as (5), (6) and 
(7) in Fig. |B]for </) = 0.38 and v s = 0. (b) Details of the unsuccessful attempt at shedding one pulse 
from the cathode for J < J c . (c) Successful nucleation of the new pulse after time (4). 
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FIG. 10. (a) Complex current signal J(t) for <p = 0.5 and a time interval At ~ 650. (b) Same 
for cf) = 0.52 and At « 400. (c) Three-pulse electric field profile for cj> = 0.52 and At 400. In all 
cases, v s = 0. In (a) and (b), the horizontal line marks the critical current J c . 
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FIG. 11. Poincare diagram depicting current maxima vs bias illustrating loss of periodic oscil- 
lations in narrow bias intervals. Parameters are r c = 10, r a = 50 and v s = 0. 
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FIG. 12. Complex current signals for v s = 0.1, r c = 10, r a = 50 and several bias values: (a) 
4> = 0.41 (periodic signal), (b) 4> = 0.411 (periodic signal with a longer period due to current 
bursts), and (c) 4> = 0.42 (aperiodic signal). In all cases, J(t) oscillates about J c = 5.77 with an 
approximate amplitude of 0.5. 
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FIG. 13. Phase plane (r w ,J) showing the nullcline dJ/dr w = (dash line) and the turning 
point for v s = and a bias <f> = 0.6. The thick line represents the trajectory of the solution for 
initial data J(0) =4.1 and r w = 26, until J c « 5.4 is reached. 
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FIG. 14. Comparison of the asymptotic expressions for J, E + and Rf, with the results of direct 
numerical solution starting with J (to = 30) = J m in ~ 2.4. Parameter values are r c = 10, r a = 90, 
4> = 0.6 and v s = 0. (a) J(t) (the horizontal line marks the critical current J c ). (b) E + (t). (c) 
Rb(t). (d) Maximum radius as a function of bias. The numerically calculated maximum radius for 
the back of the pulse is compared to r w T and r wc . The difference of the latter with the numerically 
calculated Rf, is the pulse width. 
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FIG. 15. Block matrix formulation of our numerical scheme to solve the equations for E and J. 
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